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Abstract 

This  note  is  concerned  with  the  formulation  of  a  damped  second  order  sys¬ 
tem  as  a  first  order  dynamical  system  on  a  product  space.  This  problem  comes 
from  the  desire  to  have  explicit  representations  of  the  infinitesimal  generator  of 
the  first  order  system  and,  in  particular,  of  the  domain  of  this  operator.  This 
analysis  is  motivated  by  the  need  to  find  specific  representations  for  Riccati 
operators  that  can  be  used  in  the  development  of  computational  schemes  for 
hyperbolic  control  problems.  The  approach  we  take  here  is  based  on  a  natural 
factorization  of  the  differential  operators  that  define  the  second  order  model. 
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1  Introduction  and  Motivation 


In  [2]-[4],  Burns  and  King  consider  feedback  control  problems  for  damped  hyperbolic 
systems.  Specifically,  they  are  concerned  with  obtaining  integral  representations  of 
the  feedback  control  law  for  purposes  of  designing  reduced  order  controllers  and  sen¬ 
sors/actuators.  The  kernels  of  these  integral  representations  are  called  functional 
gains.  In  [4],  Burns  and  King  use  information  about  the  spatial  support  of  functional 
gains  to  guide  the  construction  of  reduced  order  controllers  for  a  nonlinear  damped 
elastic  system.  In  order  to  use  this  information,  it  is  important  to  have  both  quali¬ 
tative  and  quantitative  information  about  these  feedback  gains.  The  very  existence 
of  these  kernels  is  not  always  obvious  and  requires  careful  analysis  of  the  system  and 
the  exact  form  of  the  damping  model.  Indeed,  the  damping  model  greatly  impacts 
these  gains  (see  [2]- [4]).  One  approach  to  modelling  second  order  damped  systems  is 
to  start  with  the  undamped  equation  and  then  “add”  a  damping  term  of  the  form 
jDox{t)  to  the  second  order  system.  Very  often  the  damping  operator  Dq  is  assumed 
to  have  the  form  of  a  fractional  power  of  the  structural  operator,  i.e.,  Dq  =  .4“  for 
0  <  o;  <  1.  This  approach  leads  to  formal  models  that  mimic  various  damping 
models  such  as  structural  damping  (a  =  1/2)  and  Kelvin- Voigt  damping  (a  =  1). 
However,  in  order  to  turn  this  formal  second  order  system  into  a  well-posed  dynami¬ 
cal  system  on  an  appropriate  state  space,  one  is  often  faced  with  having  to  deal  with 
fractional  powers  of  differential  operators  leading  to  pseudo-differential  operators.  In 
this  note  we  present  several  formulations  of  this  problem,  one  of  which  makes  use  of 
physics  based  modelling.  In  many  cases  this  approach  can  greatly  simplify  the  anal¬ 
ysis  and  the  resulting  first  order  system  has  an  explicit  representation  that  avoids 
pseudo-differential  operators. 

2  Abstract  Second  Order  Damped  Models 

Let  iy  be  a  Hilbert  space.  We  assume  that  .4  is  a  self-adjoint,  strictly  positive  operator 
on  H  with  domain  Dom{A)  dense  in  H.  Consider  the  undamped  second  order  control 
system 

-I-  Ax{t)  =  Bu{t),  (2.1) 

where  H  is  a  compact  linear  operator.  It  is  well  known  that  in  this  case,  (2.1)  is  not 
stabilizable  (for  example,  see  [6]).  All  elastic  systems  have  some  internal  damping  and 
the  exact  form  of  this  damping  is  important  in  the  analysis  and  solution  of  control 
problems  for  elastic  systems.  In  this  paper  we  concentrate  on  the  development  of 
explicit  state  space  models  for  the  uncontrolled  systems.  The  application  of  these 
inodels  to  control  design  will  appear  in  a  future  paper. 

The  general  mathematical  model  often  used  as  a  prototype  for  describing  con¬ 
trolled  elastic  systems  with  internal  damping  is  obtained  by  adding  a  damping  term 
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(2.2) 


of  the  form  ^A^x{t)  to  (2.1)  producing  the  abstract  equation 

x{t)  +  Ax{t)  +  ^A°‘x{t)  =  Bu{t), 


where 


0  <  7,  0  <  a  <  1. 

To  write  this  second  order  system  as  a  first  order  dynamical  system  one  first  defines 
the  spaces 

V  =  DoTn{A^^^)  (2.3) 


and 


with  inner  products 


and 


E  =  Dom{A^^^)  xH  =  VxH,  (2.4) 

(«i,  U2)v  =  (2-5) 

(^A^^'^Ui,A^^‘^U2)^+{Wi,W2)h  =  {uuU2)v  +  {'^i^'^^)h^  (2-6) 


respectively. 

Let  Aa  denote  the  operator  defined  on  E  by 


0  / 
-A  -7^“ 


with  domain  Dom{Aa)  defined  by 

Dom{Aa)  =  Dom{A)  x  ^Dom{A^^^)  n  Dom{A°‘)  . 

It  is  well  known  (see,  [1]  and  [5])  that  the  operator  A^  is  densely  defined  and 
dissipative  on  E  and  hence  closable.  The  closure  of  Aa,  denoted  by  Aa,  generates  a 
strongly  continuous  semigroup  on  E  and  the  domain,  Dom{Aa),  is  given  by 

e  E  \  w  eV,  {Au  +  'yA°‘w}  €  ifj  .  (2.7) 

Although  (2.7)  provides  one  representation  of  the  domain  of  Aa,  other  (more  explicit) 
representations  can  also  be  obtained.  We  proceed  to  describe  two  formulations  of  this 
closure  based  on  factorizations  of  A  which  we  denote  as  and  Aa,2-  In  order  to  keep 
this  note  at  a  minimal  length  and  yet  present  the  basic  ideas,  we  restrict  ourselves  to 
the  case  where  1/2  <  a  <  1.  In  this  case  the  following  result  provides  a  representation 
of  Dom{Aa)  in  terms  of  the  fractional  powers  of  A.  The  first  formulation,  Aa,i,  is  a 
special  case  of  Theorem  1.1  in  [5]. 


Dom{A 


a)={ 


2 


Theorem  2.1  If\<a<  I,  then  Dom{Aa)  -  Dom{Aa,i)  where 


Dom{Aa,i)  = 


and  if  z  e  Dom{Aa,i),  then 


e  E  \  u  e  Dom{A^^'^  “),  w  eV, 
{A^-'^u  +  jw}  €Dom{A^)}, 


Aa,lZ  A(jii 


u 

w 

w 

— A“{A^“®u  +  jw} 

Moreover,  Aa,i  generates  an  analytic  semigroup  on  E. 


(2.8) 


(2.9) 


One  way  to  view  the  representation  (2.8)-(2.9)  is  to  think  of  “factoring”  (2.2)  so 
that  it  is  written  in  the  form 

£(t)  +  A°‘{A'-~°‘x{t)  +  jx{t)}  =  Bu{t),  (2.10) 

and  then  constructing  the  first  order  system  based  on  this  model.  However,  this 
approach  does  not  always  capture  the  true  physics  of  the  problem.  In  the  next 
section  we  consider  a  second  factorization  and  compare  the  corresponding  first  order 
model  to  the  system  with  Aa,i  given  by  (2.8)-(2.9). 


3  A  Symmetric  Factorization 


Another  factorization  of  (2.2),  that  seems  equally  justified,  is  based  on  factoring  A 
as  A  =  and  writing  (2.2)  in  the  form 

x{t)  +  A^/'^{A^I‘^x{t)  +  7A"“^/^i:(t)}  =  Bu{t).  (3.1) 

This  is  a  very  natural  factorization  and  it  leads  to  a  “physics  based”  formulation  of 
the  first  order  model.  This  form  of  (2.2)  leads  us  to  consider  the  operator  Aq,2  defined 
on  E  by 


DomiAaa) 


^  E  \  u,w  £V,  +  7A®  G  V 


and  for  2  G  Dom{Aa,2), 


(3.2) 


■AafiZ  —  Aa^2 


W 

— A^/^{A^/^u  +  7A““^/^'u;} 


(3.3) 


At  first  glance  it  might  appear  that  Aa,i  and  Aa,2  have  little  in  common.  However, 
the  following  result  establishes  the  equivalence  between  Aq,i  and  Aa,2  for  the  case 
where  1/2  <  a  <  1.  Hence,  it  follows  that  (3.2)-(3.3)  provides  another  representation 
of  the  closure  of  Aq.. 
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Theorem  3.1  If  ^  <  a  <  1,  then  Aa  =  Aa,i  =  Aa,2- 

Proof.  If  ^  <  Q;  <  1,  then  0<a-5<|and5<|-O!<l.  It  follows  that 

g  Dom{A^^^)  C  Dom{A‘^-^/^).  (3.4) 

Let  z  =  [u,w]’^  e  Dom{Aa,i).  We  show  that  2  G  Dom{Aafi)  and  Aa,iz  =  Aa,2Z- 
Since  z  =  [u,wY  €  Dom{Aa,i),  it  follows  that 

u  €  Dom{A^I'^~°')  C  Dom{A^/‘^),  w  e  Dom{A^^‘^)  C  Dom{A°'~^^^) 

and 

{A^~°‘u  +  7iy}  €  Dom{A°'). 

If  y  =  -A‘^{A^~°‘u  +7u;}  6  H,  then  A-'°‘y  G  Dom{A°‘),  where 

A~^y  =  -{A^~°‘u  +  7ty}  =  -A^^'^{A^^‘^~°‘u  +  'yA~^^'^w}.  (3.5) 

Also, 

=  A-^^^-^y  =  A-^{A-^/^y)  G  Dom(A“).  (3.6) 

Since 

j^i/2-au  —  a~°'{A^^^u)  G  Dom{A°‘), 
it  follows  from  (3.6)  that  G  Dom{A°‘).  Hence, 

-{A'/2^  +  7A“-'/2^}  =  + ^  g  Dcym{A^I%  (3.7) 

and 

_^i/2{^i/2^  +  7A“-'/2«,}  =y£H.  (3.8) 

Combining  (3.4),  (3.7)  and  (3.8),  it  follows  that  z  G  Dom{Aa^2)  sad  Aq.^2^  =  A^^iz. 
Conversely,  assume  that  z  =  [u^wY  ^  Dom{AaY)-  Let 

y  =  -A^!^{A^I\  +  7A“-'/'u;}  G  H 

and  observe  that 

A-y^y  =  -{A^l-^u  +  7A“-'/"u;}  =  +  ^A-^'^w}  G  L>om(A'/2),  (3.9) 

implies 

+  7A-'/2w}  =  A-'/'-“y  =  ^-'/^(A-^y)  G  Dom(A'/2).  (3.IO) 

However,  ^A~y‘^w  G  DomiA^f^),  so  it  follows  from  (3.10)  that 

Ai/2-a^  G  Dom(A^/^). 
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Hence, 


_^l/2{^l/2  cc^  1/2^1  ^  ^  (^““y)  G  Dom{A^),  (3.11) 

and 

+  =y  G  H".  (3.12) 

Note  that  w  £V  =  Dom{A^I^)  and  {A^-^^u  +  7^;}  =  6  Dmn{A°‘).  Hence, 

A^~°‘u  =  A~°'y  -  7w  G  Dom{A^/^).  If  y  =  A^^‘^{A^~°‘u),  then 

u  =  A'^-\A-^/‘^y)  =  ^“-3/2y  =  ^-(3/2-a)^  ^  Dom(^3/2-“). 

Combining  (3.4)  with  (3.11)  and  (3.12),  it  follows  that  2  G  Dom{Aa,\)  and  that 
Aa,\z  =  AafiZ.  This  completes  the  proof. 

This  result  shows  that  there  are  several  equivalent  representations  of  the  closure 
of  .4^.  Although  (2.8)-(2.9)  and  (3.2)-(3.3)  both  provide  explicit  characterizations  of 
the  domain  of  this  closure,  both  characterizations  are  in  terms  of  fractional  powers 
of  A.  However,  the  representation  (3.2)-(3.3)  can  be  very  useful.  For  example  the 
following  theorem  follows  from  a  direct  calculation. 


Theorem  3.2  /f  |  <  o;  <  1,  then  the  Hilbert  adjoint  [^^^2]*  is  defined  on  the  domain 


Dom{[Aa,2]*) 


£  E  \  u,w  {A^/'^u  —  7^“  G  V 


(3.13) 


by 


[Act, 2]' 


[Aq,2. 


* 


u 

—w 

w 

(3.14) 


Observe  that  although  Aq,,2  and  [^0,^2]*  have  a  similar  structure,  they  have  diflferent 
domains.  In  particular,  Dom{[Aa,2]*)  #  Dom{Aa,2)  and  Dom{[Aa,2]*)^D(m{Aa,2)  C 
Dom{A)  X  Dom{A'^).  Again,  the  adjoint  is  given  in  terms  of  fractional  powers  of  A, 
and  in  many  cases,  these  operators  are  pseudo-differential  operators  without  simple 
explicit  representations.  A  third  approach  to  this  problem  is  based  on  returning  to 
fundamental  physics. 

4  A  Physics  Based  Factorization 

The  factorization  (3.1)  is  a  special  case  of  a  more  general  form.  If  we  define  the 
operators 

S  =  A^/2  and  D  = 

then  (3.1)  can  be  written  as 

x{t)  +  S*{Sx{t)  -f  'yDx{t)]  =  Bu{t).  (4.1) 
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Observe  that  S*  =  S  =  since  A  is  assumed  to  be  self-adjoint  and  positive 
definite.  Moreover,  the  basic  spaces  given  in  (2.3)  -  (2.6)  are  defined  in  terms  of 
Dom{S)  by 

V  =  Dom(A^^^)  =  Dom(S),  (4.2) 

and 

B  =  Dom(S)  xH  =  VxH,  (4.3) 

with  inner  products 

{ui,U2)v  -  {Sui,  Su2)jj  (4.4) 

and 

\  =  {SUi,SU2)fi  +  {Wi,W2)ff  =  {Ui,U2)v  +  ^  ('^•^) 

/  E 

respectively. 

However,  it  is  sometimes  more  useful  to  use  a  different  factorization.  We  illustrate 
the  basic  idea  by  restricting  attention  to  a  simple  ID  wave  equation.  Although  the 
presentation  here  is  focused  on  this  example  problem,  the  approach  can  be  extended 
to  a  wide  variety  of  2D  and  3D  problems  in  elasticity.  Consider  the  problem  of  an 
undamped  vibrating  string  on  the  interval  0  <  s  <  1  with  fixed  left  end  and  free  right 
end.  If  w{t,s)  denotes  the  displacement  of  the  string  and  a{t,s)  denotes  the  stress, 
then  the  wave  equation  becomes 

(4-6) 

with  displacement  boundary  condition  at  s  =  0 

w{t,0)  =  0,  (4.7) 

and  “natural”  boundary  condition  at  s  =  1 

a{t,  1)  =  0.  (4.8) 

The  strain  is  defined  by  £{t,s)  =  §;w{t,  s)  and  if  one  uses  the  stress-strain  law 

a{t,s)  =  Te{t,s),  (4.9) 

then  the  equation  (4.6)  becomes 

(4.10) 

The  appropriate  boundary  conditions  are 

w{t,  0)  =  0,  1)  =  0.  (4.11) 
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(4.12) 


On  the  other  hand  if  one  uses  a  dynamic  stress-strain  law  such  as 

d 

a{t,  s)  =  T€{t,  s)  +  s), 

then  the  equation  (4.6)  becomes 

d  (  d  1 

'  S  ra;'"*'- *>r 

In  this  case  the  appropriate  boundary  conditions  are  now 

w{t,  0)  =  0,  a{t,  1)  =  1)  -f  1)  I  =  0- 


52 

dt'^' 


(4.13) 


(4.14) 


The  partial  differential  equation  in  (4.13)  with  boundary  conditions  (4.14)  can  be 
written  as  a  second  order  system  in  the  space  of  virtual  displacements  H  =  L‘^{0, 1). 
We  also  introduce  the  space  of  deformations  (or  strains)  E  =  L‘^{0, 1)  and  define  the 
operator  S  from  H  into  E  on  the  domain 


Dom{S)  =  Hl{0, 1)  =  {«;(•)  e  H\0, 1)  |  u;(0)  =  o}  (4.15) 


by 


[5u;(-)](s)  = 


Then  the  adjoint  of  S  is  defined  on  E  into  H  by 

Dom{S‘)  =  Hl,(0, 1)  =  {<t(.)  €  K*(0, 1)  |  a(l)  =  o} 


(4.16) 

(4.17) 


and  for  (j(-)  E  Dom(S*) 

1SV(.)1(S)  =  -£<7(s).  (4.18) 

If  the  damping  operator  D  :  i?i(0, 1)  -)•  ^^(o,  1)  is  defined  by  D  =  S,  then  the  wave 
equation  (4.13)  with  boundary  conditions  (4.14)  can  be  written  as 


x(t)  -h  S*{rSx(t)  +  7-Di(^)}  =  0.  (4.19) 

Observe  that  we  do  not  distribute  S*  through  the  brackets.  In  fact,  (4.19)  is 
the  proper  form  of  the  physics  based  second  order  model  for  Kelvin- Voigt  damping. 
Moreover,  the  first  order  form  of  (4.19)  is  easily  expressed  in  terms  of  the  basic 
operators  S  and  D.  To  construct  the  first  order  model  we  define  the  spaces 

V  =  Dom(S)  =  Hl(0, 1)  (4.20) 


and 

B  =  V  X  N  =  Dom(S)  X  E  =  Hl(0,l)  X  (4.21) 
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with  inner  products 
and 


(wi,M2)v  =  {tSui,Su2) 


H 


(4.22) 


Ui 

^2 

.  "^1 . 

5 

-  '^2  _ 

{tSuuSu2)jj  +  {Wi,W2)u  =  {Ui,U2)v+  {'Wi,W2)h,  (4.23) 


J  /  E 


respectively.  Now  define  ^  on  £■  =  1)  x  1)  by 


u 


w 


Dom(A)  =  = 

where  for  z  G  Dom{A), 

Az  —  A 


^  E  \  u,w  E  V,  {rSu  +  'yDw}  G  Dom{S*) 


(4.24) 


U 

w 

W 

—S*{tSu  +  yJDiy} 

(4.25) 


Note  that  (4.24)  -  (4.25)  is  similar  to  (3.2)  -  (3.3)  where  a  =  1.  However,  there 
are  two  primary  differences  between  the  operators  A  and  Ai^2-  The  operator  S  is  not 
a  fractional  power  of  A  =  S*S  and  S*  ^  S.  Therefore,  the  operator  defined  by  (4.24) 
-  (4.25)  is  (in  general)  less  “complicated”  than  Ai^2  and  yet  we  still  have  the  following 
easily  established  result. 


Theorem  4.1  The  operator  A  defined  by  (4-24)  -  (4-25)  generates  an  analytic  semi¬ 
group  on  E  =  Hl{0^  1)  X  L^(0, 1). 


The  above  theorem  has  been  extended  to  other  PDE  based  models  of  elastic 
systems.  Currently,  we  are  working  on  a  framework  that  applies  to  general  abstract 
second  order  systems  of  the  form  (4.19).  This  framework  has  the  advantage  that  the 
underlying  spaces  and  operators  are  basic  differential  operators  defined  on  standard 
Sobolev  spaces.  In  addition,  the  physics  based  factorization  is  the  natural  choice 
when  developing  approximations  (see  [1]). 

It  is  interesting  to  observe  that  if  one  starts  with  the  undamped  equation  (4.10) 
with  physical  boundary  conditions  (4.11)  and  simply  “adds  a  damping  term”,  then 
it  is  possible  to  lose  the  correct  physical  boundary  conditions.  For  example,  define  A 
on  H  =  T^(0, 1)  with  the  domain 


by 


Dom{A)  =  |u;(-)  G  H^O,  1)  ;  u;(0)  =  0,  =  o|  (4.26) 


(4.27) 
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The  wave  equation  (4.6)  becomes 


x{t)  +  Ax{t)  =  0,  (4.28) 

and  if  one  adds  a  damping  term  with  Dq  =  (i.e.  a  =  1),  one  obtains  the  second 

order  system 

x{t)  +  jAxit)  +  Ax{t)  =  0.  (4.29) 

Note  that  (4.29)  is  actually  the  abstract  form  of  the  damped  wave  equation 

+  ■y—w{t,  *)  I  =  0,  (4.30) 

with  boundary  conditions 


tt;(t,0)  =  0,  T—w{t,l)  =  0. 

However,  the  correct  physical  boundary  conditions  (given  by  (4.7)  and  (4.8))  should 
be 

w{t,  0)  =  0,  1)  +  1)  I  =  0.  (4.31) 

Therefore,  the  system  (4.29)  is  not  the  abstract  form  of  the  physics  based  model 
defined  by  the  partial  differential  equation  (4.30)  with  boundary  conditions  (4.31). 
We  close  by  remarking  that  although  (4.29)  does  not  “capture”  the  correct  physical 
boundary  conditions,  that  is  not  to  say  that  Aa,i  defined  by  (2.8)-  (2.9)  is  not  im¬ 
portant  in  the  study  of  such  systems.  However,  it  is  crucial  to  understand  that  this 
system  may  not  be  the  abstract  form  of  the  physical  problem  that  is  under  control. 

5  Conclusions 

In  this  paper,  we  present  three  formulations  of  the  abstract  form  of  damped  second 
order  systems  based  upon  different  factorizations  of  the  structural  operator.  One  form 
which  is  based  upon  the  physics  is  especially  useful  in  that  the  underlying  operators 
are  differential  operators  with  simple  explicit  representations.  Further,  this  form 
captures  the  physics,  specifically,  the  correct  boundary  conditions.  This  formulation 
has  been  extended  to  other  PDF  based  models  of  elastic  systems.  Additionally,  we 
are  currently  working  on  a  framework  that  applies  to  general  abstract  second  order 
systems  of  the  form  (4.19).  This  framework  has  the  advantage  that  the  underlying 
spaces  and  operators  are  basic  differential  operators  defined  on  standard  Sobolev 
spaces.  Moreover,  this  framework  is  a  natural  choice  when  developing  approximations 
(see  [1]). 
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